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Abstract 

Recently, we have reported theoretical studies (J. Chem. Phys. 129, 054703, 2008 and J. Chem. 
Phys. 130, 086101, 2009) on the rate of energy transfer from an electronically excited molecule to 
graphene. It was found that graphene is a very efficient quencher of the electronically excited states 
and that the rate tx (distance) -4 . The process was found to be effective up to 30 nm which is well 
beyond the traditional FRET limit. In this report, we study the transfer of an amount of energy 
Ml from a dye molecule to doped graphene. We find a crossover of the distance dependence of the 
rate from (distance) -4 to exponential as the Fermi level is increasingly shifted into the conduction 
band, with the crossover occurring at a shift of the Fermi level by an amount Ml/2. 
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I. INTRODUCTION 



Excitation energy transfer involving carbon based materials is interesting due to the fact 
that using such materials, it is possible to measure distances well beyond the traditional 
FRET limit. In our earlier papers, we have analyzed the process of resonance energy transfer 
from an excited dye molecule to a sheet of graphene lr|3|- The rate was evaluated as a 
function of the distance z of the molecule from the graphene sheet. We have found the 
process of energy transfer to be very efficient and the rate has a z~ A dependence on the 
distance. Our report was the first study on energy transfer to graphene. Recent experiments 
that have been performed after our theoretical studies have infact found efficient energy 
transfer to graphene and the process was found to be useful in identifying graphene flakes 
both on substrates and in solution |4|. Quenching by graphene was also found to be useful in 
obtaining good resonance Raman signals from fluorescent samples (5)], in fabricating graphene 
based devices 6] and in quantitative DNA analysis [7, 8]. We have also studied the process 
of energy transfer from fluorophores to carbon nanotubes and found a d~ 5 dependence 
Quantum chemical studies on energy transfer involving two carbon nanotubes have also been 
reported 10|. All the above studies involve energy transfer to extended charge densities 
of carbon based materials and hence found a deviation from the (distance)' 6 dependence, 
which has been obtained within the dipolar approximation. Such deviations from the dipolar 
approximation have also been found in polymers 111 ], quantum wells 12] etc. The Fermi 
surface of undoped graphene is a set of six points known as the K-points. As a result of 
this, the density of states at the Fermi level is zero. It is possible to shift the Fermi level 
of g raphene away from the K-point experimentally, either by electrical or chemical doping 
13^ll4|. This will make the density of states at the new Fermi level non-zero. In this letter we 



study the effect of shifting the Fermi level on the distance dependence of the rate of energy 
transfer to graphene. We imagine that the Fermi level is shifted into the conduction band 
to a level with magnitude of wave vector, hp. To keep the calculations simple, we use the 
Dirac cone approximation, which allows us to get analytical expressions for the rate at large 
distances. We note that as we are shifting the Fermi level by rather large amounts, there 
will be sizeable corrections to the rate due to deviations from the Dirac cone approximation 
and hence our conclusions are of qualitative nature. 
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II. MODEL FOR THE RATE 



We consider the process of excitation energy transfer from a dye molecule to doped 
graphene. Since the energy donor (dye molecule) has a localized electronic charge density, 
we think of the interaction between the donor and the acceptor as that between the transition 
dipole of the donor, zz^, given by fx® g = — e J dri?p®* (ri) (ri) and the transition charge 
density p (r 2 ) = —eip^* (r 2 ) tpt ( r 2) of the acceptor [3|, |9j. The matrix element for interaction 
is given by 

U = ii%- V$, (1) 

where $ is the electrostatic potential at the point r (the position of the donor) due to the 
charge density p (r 2 ) and is given by 

* {I) = J- [ dt2 ^M- . (2) 
Aire J |r — r 2 | 

As a result of energy transfer, an electron in graphene with wave vector kj is excited to a 
level with wave vector ky. We define kj = k, + q, where qh is the momentum transferred 
to graphene. The rate of energy transfer can be evaluated using the Fermi golden rule and 
is given by 

^ = ¥ EE i u *« r kk* - ^ - ( 3 ) 

k 4 q 

We use the tight binding wave functions of graphene and evaluate the matrix element (for 
details, see [l 3 ] ) using Eqs. ([I]) and (J2J) and find it to be 

U = 4^4 [ ei(5kt+q ~ 5kl) ~ !] A*S ' + k) e- qz e~ icl - x , (4) 

where q = - is the unit vector in the direction of q and k is the unit vector in the z direction. 
We have also used r = (X, z), with X being parallel to the plane of graphene. A is the area 
of the graphene lattice and 

4 = tan -1 I -M = ip k , (5) 



where (f-^ is the angle that the vector k makes with the x-axis. We substitute Eq. (]1J) for 
the interaction energy into Eq. ([3]) to get 

Ahe^A 2 



k = JKh~2 E K-^ + fc ) I' exp(-2gz)G(q), (6) 



q 
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where 

G (q) = J2 [1 - cos foc-m - V?kJ] 5« +q - ^ - tiSl). (7) 

ki 

When the Fermi level of graphene is shifted into the conduction band, the rate of energy 
transfer has contributions from two different sets of transitions. In the first, kj lies in the 
valence band with < ki < oo and kj lies in the conduction band with kp < kf < oo. In 
the second, both kj and kj lie in the conduction band with < ki < kp and kp < kf < oo 
(see Fig. [I]). The total rate can thus be written as a sum total of both the contributions, 
k — k\ + ki. k\ and ki are both given by 

+ k)| 2 exp(-2g^)G l (q), (8) 



ire 



Ahe 2 A 2 



but with differing expressions for Gi (q). 

G 1 (q)= J2 [I ~ cos (tp^-tp^SiE+^-E^-m^dki + dl-kp). (9) 

k;S valence band 

The theta function is introduced to satisfy the condition, kp < kf < oo. In a similar 
fashion, 



Gi (q) = [1 + cos (^ ki+q - </? k J] 6{E+^-E+-m)Q (\ki + q| - kp) 6 (kp - |k t 

k;6 conduction band 

(10) 

The two theta functions satisfy the conditions, < ki < kp and fcp < kf < oo. We now 
evaluate C?i (q) and G2 (q) separately. We replace the sum over kj in the expression for 
G\ (q) by an integral and use the linear dispersion relation for the energy levels of graphene 
(E± = ±v f k) to get 



Gi (q) 



.4 



4:7r 2 Vf 



dki 



1 - 



(kj + q) 



(ki + q) 



5(\ki\ + |kj + q| )0 (|ki + q| - kp) . (11) 

v f 



Introducing a new variable kj defined by kj = kj 



Gi (q) 



1 - 



V— 3 

I ! 2 I 



I leads to 



* ( ik;- f 1 



|kl + f| 



x 



e(|k; + || -kp). 



(12) 



We choose the direction of q as the x-axis and then make another change of variable to r 



given by r 



^1 to s et 



Gi (q) 



Aq 2 



16lT 2 Vf 



JdL 



1 - 



e 



5 


U 


r — i 


+ 


r + i 


\ _ m 










) v f . 



X 



r + 1 
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When Femii level is shifted 
into the conduction band 



FIG. 1: A schematic of the energy bands, showing the Fermi wave vectors, the initial and the final 
wave vectors corresponding to excitation energy transfer. Ep = vjkp is the location of the new 
Fermi level on doping graphene. 



Using r = (x, y), the above equation can be rewritten as 



oo oo 



Gi(q) = i^ / dxjdy 



-oo 



1 - 



x 2 +y 2 -l 



( X _iy + J ,2 + n x + l y + y 2 _^ 







X 



2 -ki 



(14) 



We now change over to elliptic coordinates defined by x = fiv and y = \V (/i 2 — 1) (1 



The transformation gives dxdy 



vV-i)(i-« 2 ) 



dfidv. With the above transformation, we 



get 



A( f f , f J ll-v 2 .f Ml\^rq, , , 
G i fa) = TJH. I "A* / dv \\— — 7° [ Qt* ~ —J e Ca* + v ) - k F 



An 2 Vf J J y /i 2 — 1 
1 -1 v 

The integral over \i can be easily performed to get 

G\ ( q ) = ^ e{m ' qVf) jdvVT^e 



(15) 



q ( Hfl 



2 V^/ 



(16) 



The above equation can be rewritten as 



Gl (q) = V e ( «i- TO) e 



47T 2 



^(ftfi) 2 - gV f 



2£ F - 

guy 



dvVT 



;i7) 



Mai[-1, 



We now substitute back the above expression into the rate expression of Eq. (jSJ) and convert 
the sum over q to an integral to get 



fci 



y/(Kl) 2 -qW f 



2E F -hQ 
qv f 



J dvy/l — v 2 x 



Max\-1, 



2tt 

jde 



(18) 



/*g • ( z'q + k 



where (g, 0) are the polar coordinates of q. After performing the integral over 9, we average 
over all possible orientations of the donor transition dipole (see {2, 3]) to get 



9 9 1 - XJ 



1 _ 2E F -nn 



f dvyl — v 2 . 



Max[-l, 2EF ' hn ] 



(19) 



Evaluation of the integral over v leads to 



96ft€ 2 7r 2 



f daa 3 e- 2<!Z e( - m - qv ^ Q 

Jmf^v* 



2E F -nn 

qv f 



[f - - « 2 + sin" 1 it)] , (20) 



where it is defined by it = Max[— 1 



We now evaluate G 2 (q), defined by Eq. ( TTOT) . Using the same procedure as before, we 
find that 



G2 (q) 



Aq 2 
4ir 2 Vf 



00 1 



dfi / rff 



2 - 1 



5 ( qv 



v f ) 



e 



9 1 r 9 

- (/i + l>) - fc F J © y$F 2 ^ 



(21) 



The integral over i> can be easily performed to get 



G2 (q) 



Aq 2 e(qv f -htt) 



An 2 



q 2 v 2 - {my { 



00 

J dfi^n 2 - le 



g / Ml 



2 Vg^/ 



/U - fci 







k F --[n 



- 



qv f J 
(22) 



The above equation can be rewritten as 



2E F +h(l 



G2 (q) 



Aq 2 Q(qv f -m) 
47f2 !q 2 v 2 - (Ml)'' 



djiy fi 2 — 10 



Max[l, 2EF ~ m ] 



2£ F + htt 

qv f 



- 1 



(23) 
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Substituting the above expression back into the rate expression of Eq. ([8]) gives 



o 



2E F +nn 



2E F +htl 

/ dlXyJ /i 2 — lX 



Max\l, 



(24) 



2vr 




D 



iq + k 



The integral can be performed easily, followed by an averaging over all possible orientations 
of the donor transition dipole to get 



2 2 00 

. _ e r j nn s r -2az e(gw/-fin) ^ 
fc 2 - 48rt£ 2, 2 J dgg e ^_^u 

The integral over can be evaluated to get 



2E F +nn 



J d[i\J /J? — 1. 



(25) 



Maxfl, 2gf, ~ ,in ] 



X 



^ q *v)-{mf 
\Jr 2 - 1 + sVs 2 - 1 + log [r + V^ 2 - lj - log [s + Vs 2 - lj ) 



(26) 



where r = Max[l, 2g ^~ r ' n ] and s = 2E g^~ f hQ ■ The integrals over g in the expressions for k\ 
and k 2 can be performed numerically, thus getting the total rate of transfer when the Fermi 
level in graphene is shifted into the conduction band. 



III. LARGE z BEHAVIOR OF ki AND k 2 

It is easy to analyze the large z behavior of k\ and k 2 (see the Appendix for the detailed 
analysis). When E F < ^ and z > with Ae = ^ — E F and Ml » Ae, the long-range 
behavior of k\ is given by Eq. ( 1321) of the Appendix as 

e 2 u 2 

1 2567iQh 2 e 2 z i ' 1 } 

In the case when Ep > ^ and z > with KVL >> Ae, 

e 2 H 2 |Ae| 3/2 , s 

fci = — —rr, e u / , 28 

96\/2fin 2 e 2 7r 3 / 2 ^ /2 z 5 /2 

as obtained in Eq. ( )39|) of the Appendix. The major contribution from the k\ term to the 
rate comes only when E F < 4p and it has a power law dependence {z~ A ) on the distance. 
When Ep > A?, the contribution from k\ decreases exponentially with z and hence is very 




FIG. 2: The rate of energy transfer as a function of distance, as the Fermi level is shifted into the 
conduction band. The emission energy of the fluorophore is taken to be Ml = 2.0 eV. 

small. The long-range behavior of &2 for both Ep < and Ep > 4p is given by Eqs. (jH|) 
and (l49l) of the Appendix as 

j — 2zfifi 

fe = f$fj;'\Zz lr*r=T - log [r + V^T] ) , (29) 

where r = 1 + 

Therefore, the large z behavior of k-z is exponential. Thus, in the case when Ep < 4p, 
the rate of transfer to doped graphene has a power law dependence on the distance (arising 
from the k\ term), while when Ep > the rate has an exponential dependence (arising 
due to both k\ and k-i terms). Therefore, as the Fermi level is increasingly shifted into the 
conduction band, there is a crossover of the distance dependence of the rate from power law 
to exponential and the crossover occurs over a region of Ep centred at 
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FIG. 3: The rate of energy transfer as a function of distance, as the Fermi level is shifted into the 
conduction band. The emission energy of the fluorophore is taken to be Ml = 2.0 eV. 

IV. RESULTS 

We performed numerical calculations for evaluating the rates of energy transfer from a 
fluorophore to doped graphene. We take the emission energy of the fluorophore, Ml to 



be 2.0 eV. Such low energy emission has been found in squarylium dyes [15j. We take 
fi eg = 4.5 D. Fig. [2] shows a plot of the logarithm of the rate as a function of the logarithm 
of the distance as the Fermi level is shifted increasingly into the conduction band in the range 
0.2 — 2.0 eV. The log-log plot is linear showing that the rate has a power law dependence 
{z~ 4 ) on the distance for Ep < 1.0 eV. After around 1.0 eV, there is deviation from linearity 
in the log-log plot and the dependence becomes exponential. In order to look into the cross 
over region more closely, in Fig. [3], we show a plot of the rate as the Fermi level is moved 
into the conduction band in the range 0.95 — 1.05 eV. This clearly shows that there is a 
crossover of the distance dependence of the rate from power law to exponential as the Fermi 
level is increasingly shifted into the conduction band. It should be possible to observe this 
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FIG. 4: The rate of energy transfer as a function of the Fermi energy of graphene, at a distance of 
10 A. The figure shows plots of k\ and k 2 separately along with the total rate of transfer. 

effect experimentally. 

Fig. H] shows a plot of the k\ and k 2 terms, and the total rate of transfer as a function 
of the Fermi energy of graphene at a fixed distance, z = 10 A. The contribution from the 
k\ term decreases as the Fermi level is increasingly shifted into the conduction band, while 
that from the k 2 term increases. From the figure, it is clear that the total rate is governed 
mainly by the k\ term up to the crossover point. Beyond that, the contribution from ki 
term is small and the total rate is governed by the k 2 term. This is clearly a density of states 
effect. 

V. CONCLUSIONS 

We have studied the process of energy transfer from a fluorophore to doped graphene. We 
have analyzed the distance dependence of the rate of transfer as the Fermi level of graphene 
is shifted away from the K-point into the conduction band. We find a crossover of the 
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dependence from power law (z -4 ) to exponential as the Fermi level is increasingly moved 
into the conduction band. The point of crossover is at a shift of the Fermi level by Ml/ 2. 
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VI. APPENDIX 



A. Behavior of rate constant, k\ 



We first look at k\. We consider two separate cases. 

1. Case I: We consider the case where E F is close to but less than it. We put 
Ae = ^ — Ep. k\ is then given by 



48fi<E 2 7T 



1+ 2Ae 

qv f 



f dvy/l — V 2 . 



Max[-l-^\ 

1 ' qv £ J 



(30) 



As Ae > 0, 6 



1+ 2A£ 

qv f 



1. For 2 > 2A - , the major contribution to the above integral 



is from q £ (0, — ). In this range 



7T 



(31) 



Afoa[-l ,-^1 



If >> Ae, then, in this range, (M7) 2 — q 2 v 2 ~ (h£l) 2 and the integral in Eq. ( 130]) 
may be approximated as 



e 2 /x 2 



967r/i 2 e 2 fi 



dqq e 



Z p -2qz 



eV 2 



2567rfifr 2 e 2 z 4 ' 



(32) 



2. Case II: We now consider the case E F > ^. Hence, ^ — Ej 



|Ae|. Then, 



48/le 2 ' 



■AM 

qv f 



f dvVl 



V* 



Ma,[-l,^f] 



For large z (^z >> 4^yj; this may be approximated as 



e 2 /4 
48fre 2 vr s 



2|Ae| 
v f 



dvyl 



v 2 . 



(33) 



(34) 



For z such that 4 l Ae l* >> 1, in the above integral over g, the contribution is from q in 



<7 



the vicinity of 3^1. We now change the variable of integration from q to y defined by 



q = 2 ^ y - Using the above transformation, fci is given by 



e 2 /i 2 3 |Aer 
3he 2 ir 2 vj 



L...,.^ e(«n-2|A e |,) (35) 

J ■ l{mf -4|Ae|V { 
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Now, the major contribution to the above integral comes from values of y ~ 1. If y is 
close to unity, 

3/2 



Hence 



|wr^^(^) . (36) 
i 

V 



tl ,?W/^ (s _ lfe ^JEdW|L. (37) 

/ ^ v(^) - 4 l Ae l y 



For ^ >> |Ae| and y in the vicinity of 1, Q(Ml — 2 |Ae| = 1 and 
(m) - 4 | Aef y 2 ~ HQ. The integral over y can now be evaluated to get: 

e 2 /i 2 |Ae| 2 - 2 i Ae ^ /2|AeU\ 



48\/2fi^ 2 e 2 7r 2 w 2 ^ 2 V v f 



For large values of z, the asymptotic form of the Bessel function K 2 {^^~^ is gi 
by # 2 (^) ~ y^e"^. Therefore, we get 



lven 



e 2 /i 2 |Ae| 3/2 ^IMi t 

~ 9QV2m 2 eH 3 / 2 v 3 f /2 z 5/2 



B. Behavior of rate constant, k 2 

We now analyze k 2 . As before, we consider two separate cases. 

1. Case I: We consider the case Ep < Therefore, 2E F -hfl<0 and hence k 2 is given 
by 

2E F + hSl 

2 2 OO 9^7 

Using the two theta functions in the above expression, k 2 can be written as 

2E F +nn 2E F +nn 



2E F +hU 
qv f 



(40) 
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We now make a change of the variable of integration from q to x defined by q = — + x. 



Using this transformation, k 2 can be written as 



fc 2 



2, ,2 



-2znn v f 



HQ 



2E F +nn 

-2zx xvf-\-h£l 



48/^ C J UX ./ X 2 v 2 +2mv x J 

V J J 1 



J dfiy/fl 2 — 1. 



(42) 



For large values of z, in the above integral over x, because of the presence of e~ 2zx 
term, only small values of x are important. Therefore, for small values of x, the above 
expression can be simplified to get 



-2znn 1+- 



^ 48V2 I ~ 1 J ■ 

1 

The integral over \i can now be performed, and the upper limit in the integral over x 
can be extended to oo to get the following expression for k 2 . 



(43) 



k, 



where r = 1 + 2 ^-. 



r W (^v^T- log [r + v^T]) 



(44) 



2. Case II: We now consider the case Ep > Therefore, 2Ep — hVt > and hence k 2 



is given by 



ft 2 - 



3„-2qz @(qvf-tin) 



2E F +hfl 
1 v f 



Max[l, 



2E F +nn _ j 

qv f 



(45) 



The above equation can be simplified to get 



h = ^ J dq^ 



3 p -2qz 



1 v f 



(46) 



We now use the same procedure as was used for evaluating the integrals in Case I to 
get 



48^ 



-2zhSl H — 



/ dfly/fl 2 - 1 / cfe^J . 



(47) 



On evaluating the integral over /x, we get 



96^ 



(— s\/ s 2 — 1 + r-y/^ 2 — 1 + log [s + a/s 2 — lj — log [r + \/r 2 — lj ) 



x 



/ dx'- 



(48) 
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where s = Max[l, - 1] and r = 1 + ^§f . For E F close to ^, s = 1. Using this and 
extending the upper limit of the integral over itooo and then evaluating the integral, 
we get 

I — 2znn 

h = ./ -ar \Z, (»V^n: - l°g [r + Vr^l] ) ■ (49) 



"y 
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